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Abstract: 

In this paper, a numerical solution for inverse 
problem involving Irregular boundary value problem using 
B -cubic spline has been presented. 

Also identification of the domain of definition is 
studied. 

I. Introduction 

The inverse problems are mathematical problems 
appears in many applications of mathematics such that 
some information about the solution of the problem is given 
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and some parameters (or conditions functions ) are unknown 
in this problem . 

The inverse problem is used in various area, such 
dynamics , engineering, statistics, numerical analysis, 
mathematical modeling, etc [7] , [1] . In this paper ,the B- 
cubic spline base has been used in approximation . 


^Preliminaries 


2.1 B - Cubic Spline Base 


We first partition the interval [-2,2] by choosing the 
length of partition equal to 1, [6] . This produce the equally 
space nodes x f =-2 + i for each i = 01,2,3,4 . 


Define B : [ - 2,2 ] >[0,l] as B-cubic spline function on 

the interval [-2, 2] satisfy the smoothness properties i.e. 

Be C 2 (-00,00) suc h that : 


(2 - x) 3 - 4(1 - x) 3 - 6 + 4(1 + x) 
(2 - x) 3 - 4(1 - x) 3 - 6x 3 
(2 - x) 3 - 4(1 - x) 3 
(2 - x) 3 


- 2 < x < -1 
- 1 < x < 0 

0 < x < 1 

1 < x < 2 


Also the B Cubic Spline Base 

B = {b 0 (x). B j (x), B 2 (x), B n (x) } , [5] in the interval [a , b ] 

which has the length of partition/, = b ~ a with nodes of 


partition x t 
following : 

B,(x) =B| 


: a 


+ i*h ) for i=0,l,2,3,...,n have been define as 


X - X- 


for i=0,l,2,...,n , [5]. 


2.2 The method of solution 

Consider the Boundary value problem (B.V.P.) 
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L u = f(x, y) i n R 

( 2.1 ) 

u = 0 on d R 

Such that 


R = {(x, y) : a < x < b,0 < y < h(x) } 

See figure 1, problem 2.1 has unique solution if h(x) is 
given . Suppose that this problem has solution w*(x,y) and 
the function h(x) is unknown. 

Now to estimate h(x), it will be approximated by linear 
combination of Cubic spline functions base 

B 0 (x), B j (x), B 2 (x), B n (x) which are linearly independent 

and smoothness . 

Let h'(x) = X a k B k (x) 

k = 1 

In this paper, the least squares method is used to 
estimate the parameters a p a 2 ,,,,,a n from 

( -> . y 

J(a , , a 2 a n ) = £ ^ \u ( x . , y . , a ) - u * ( x jt y ) 

i j V 

Where 

u represents the given exact solution, 

u represents the Numerical solution of 

problem. 

— ^ 

d = (a j , a 2 , , , , , a n ) the vector of the unknown parameters 

The values of parameters can be found by minimizing J 
with respect to a, for i=l,2,...,n 

The method solution of inverse problem will be 

transform into the unconstrained nonlinear programming 

— > 

problem min J(a) 

There are many methods to solve the problem ( 2.1 ), 
the variational method has been used, [6] also the B-cubic 


\ xr 
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spline base has been used to approximate the solution of 
boundary value problem ( 2.1 ) . 

2.1.1 The Variational Method 

Let the boundary value problem : 

v 2 u =f in R 

u =0 on 3R 

Such that 


R = {(v, y):a<x<b,c<y< h(x) } 

The equivalent functional to this problem is : 

b h{x) 

J(u) = J \{u 2 x + u 2 y + 2 uf )dxdy ( 2.2 ) 

a 0 

j _ b - a 

[6] . Define the norm of x-interval as follow : ' l = 

n 

_ h(x) 

And the norm of y-interval is : K x ~ such that 

m 

n and m are two positive integer numbers [4] . 

The B cubic spline base in x has been used such that 
satisfy the boundary conditions of problem . It means that : 

B, (a, h) = B ; ( b , h) = 0 for i=0,l,2,..,m 

Also, the B -cubic spline base in y has been used such 
that satisfy the boundary conditions of problem . It means 
that : 

B;(/i(x), k x ) = B ; (0, k x ) = 0 for i=0,l,2,. . .,n 
So that . the solution is 

N,M 

U (x, y ) = ^ a tj B t (x, h)B j (y, k x ) 

i,j = 0 

Which Satisfy B.C. of problem [3]. 


(2.3) 
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Substitute equation ( 2.3 )in the equation ( 2.2 ) and 
use partial derivative with respect to parameters that 
mean : 

-^-=0 for k = 0 , 1 , 2 ,..., n ,1 = 0,1,2,..., m 

dz kl 

The Guassin method is used to solve the system of 
linear equations Az=b . 

To solve the unconstrained nonlinear programming 
.The Hooke and Jeeves method is used, Results are given in 
figure (2) which illustrate the results in case 
h(x)=l+x(x-l) 
u(x,y)=xy(l-x)(l+x(l-x)-y) 
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